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basis of equilibrium statistics. In particular, the influence of a flat hard wall on the 
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dependent structure function in position space are determined from an expansion 
in terms of the eigenfunctions of the one-particle Hamiltonian. For the bulk of the 
system, the intermediate scattering function and the dynamical structure factor 
are derived by taking successive Fourier transforms. In the vicinity of the wall 
the time-dependent coherent structure function is found to decay faster than in 
the bulk. For coinciding positions near the wall the form of the structure function 
turns out to be independent of the orientation of the wall. Numerical results are 
shown to corroborate these findings. 



PACS numbers: 05.30.Fk,75.20.-g 



Submitted to: J. Phys. A: Math. Gen. 



Time correlations in a confined magnetized free-electron gas 



2 



1. Introduction 

The equilibrium properties of magnetized quantum electron gases are partly 
determined by physical processes near its boundaries. For instance, the Landau 
diamagnetic effect of a magnetized free-electron gas is a result of currents flowing 
near the edge Hence, it is important to study the influence of a wall on both the 
static and the dynamic behavior of these systems. In a recent paper |2] we investigated 
the static correlations in a fully degenerate magnetized free-electron gas in the vicinity 
of a wall. In the bulk the decay of the pair correlation function for large separations 
depends on the orientation of the position difference: for directions orthogonal to the 
field it is exponentially fast, whereas it is algebraic for directions parallel to the field. If 
a hard wall parallel to the field is present, the picture changes: the asymptotic behavior 
of the pair correlation function for transverse directions changes from exponential to 
algebraic for any two positions close to the wall. 

Further insight in the influence of a wall on the correlation properties is obtained 
by studying dynamical correlations, which are the subject of the present paper. In 
the bulk these time correlations yield, after suitable Fourier transforms, the dynamic 
structure factor, which has been studied quite extensively for various systems, both 
theoretically and experimentally. For an unmagnetized degenerate free-electron gas 
its properties were determined several decades ago 0. The generalization to finite 
temperatures is available as well 0]. The dynamic structure factor for a magnetized 
degenerate free-electron gas is somewhat more complicated, as it involves a sum over 
Landau levels. If a wall is present as well, the translation invariance of the system is 
lost, so that taking a Fourier transform is inconvenient. Instead, it is more suitable to 
analyze the dynamical correlations by means of the structure function in configuration 
space. We will focus on its time dependence, and, in particular, on its decay for large 
time. As a starting point we shall employ the expansion of the space-time-dependent 
structure function in terms of the single-particle eigenfunctions. 

The paper is organized as follows. In the next section we shall define the structure 
function and the associated time correlation functions. In section 3 we shall derive 
expressions for the bulk space-time-dependent structure function of a magnetized 
degenerate free-electron gas. Some numerical results for this function will be presented 
as well. The sections 4 and 5 contain an analysis of the intermediate scattering 
ftmctions and the dynamical structure factor for the bulk, again with some numerical 
illustrations. In section 6 we shall verify that in the limit of vanishing magnetic field 
our results reduce to the well-known expressions for the field-free degenerate free- 
electron gas. Edge effects are the subject of sections 7 and 8. Both the cases of a wall 
perpendicular to and parallel with the magnetic field will be discussed. The analytical 
results will be checked numerically. 

2. General properties of quantum time correlation functions 

The equilibrium time-dependent structure function in position space is defined as 

5(r,r',i) - (e'«*^t(r)^(r)e-'^* i^H^'Hir')) - (^t(r)V;(r)) (V-^rOV^r')) (1) 

with tp(r) the field annihilation operator and H the many-body Hamiltonian. The 
brackets denote an equilibrium average in the grand-canonical ensemble, with the 
inverse temperature /3 and the chemical potential fi. We have chosen units such that 
h drops out. 
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For independent particles with Fermi statistics the structure function can be 
expressed in terms of single-particle wavefunctions as 

S{y, r', t)^YlY^ efi(E,^t.) _,. 1 ^ ^ QP(E^,~t,) , ^ 

3 j' 

x^*(r)^,(r')^,'(r)¥.*,(r')e'(^^-^.')' (2) 

where ipj (r) are the eigenfunctions of the single-particle Hamiltonian with eigenvalues 
Ej . The structure function can be written as the sum of an incoherent and a coherent 
part IS]- These are given as 

S,nc{r,r',t)^G;^T{r,r',t)G{r,r',t) (3) 

Scoh{r,r',t)^~\G,Mr,r',t)f (4) 

Here we introduced the time correlation functions G^^t and G, which are defined as 

G,Mry,t) . ^-^^^J-^ ^^.(r)^*(r') e"*^^* (5) 
j 

G(r,r',t)=^^,(rV*(r')e-'^^* (6) 

j 

The first function depends on fi and T, whereas the second does not. In the static 
case of vanishing t the correlation function G reduces to 6{r — r'). 

Our main interest in the following will be in the fully degenerate case with T ~ 0. 
In that case G^^t reduces to 

G,,T=o{r, r', t)=Y,9{^,~ E,) ^,(r)^*(r') e"'^^* (7) 
J 

with 9 the step function. The time correlation function for finite T can easily be 
recovered from that of the fully degenerate case, since once has 

G,.T(r, r', t) ^ d^' ^ ^ ^ G,,^=o(r, r', t) (8) 

The subscript T = will often be omitted in the following, so that G^_t=o will be 
written as G^. 

In uniform systems, or in the bulk of a confined system, the structure function 
^(r, r', t) depends on the position difference r — r' only. In this case the intermediate 
scattering function follows by taking the spatial Fourier transform of S: 

5(k, t)^ j d{r - r') e-"' ^'"^"^') S{r, r', t) (9) 

and similarly for the incoherent and coherent parts. A further Fourier transform leads 
to the dynamic structure factor: 



S{k,Lu) = J dte"^' S{k,t) (10) 
Substituting jSJ and © in @ and performing the Fourier transforms one gets 



e/3(iSj-A') + 1 

j j' 



X / d(r - r') e-^-C-') ^Kr)^, (r')^,' (r)^* (r') (11) 
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A similar expression can be written down for the coherent part. Comparing the two 
expressions one easily finds the relation 

Scoh{K t^) = ^ ^mc(k, uj) + _ ^ S'™c(-k, ~UJ) (12) 

This relation also follows from the principle of detailed balance [H] , if the symmetry 
of the coherent dynamic structure factor under a change of sign of both k and uj is 
taken into account. 



3. Time correlations in the bulk 



We consider a system of charged particles which move in a uniform magnetic field 
directed along the z-axis. The interaction between the particles is neglected. To 
describe the magnetic field we adopt the Landau gauge, with vector potential A = 
{0,Bx,Q). The single-particle Hamiltonian reads 

H = -^A + iBx-^ + ^B'^x^ (13) 

Units have been chosen such that the charge and the mass of the particles drop out, 
while c (andS) have been put to 1 as well. From now on we will often measure distances 
in terms of the cyclotron radius l/\/S. To that end we introduce the reduced variables 
^ = \^x, rj = \fBy and C — VBz. 

The normalized eigenfunctions of the Hamiltonian are: 

rjl/4 

We imposed periodic boundary conditions in the y- and z-direction, with a periodicity 
length L. The wavevector components are given by ki = 27rni/L, with integer 
for i = y, z; the reduced wavevector components are Ki = ki/VB. The Hermite 
polynomials 7J„ carry the non-negative integer label n. The associated eigenvalue is 
the sum of the kinetic energy for the z-direction and the Landau-level energy for the 
transverse directions: 

En^^n=B{n + ^) + ^kl (15) 

The time correlation function for the fully degenerate system follows by 
inserting the eigenfunctions and eigenvalues in Q. For large L the summations over 
Hi may be replaced by integrations over k^. As a result we find: 

n3/2 °o 1 
n=0 

We introduced the reduced time r ~ Bt, and the reduced chemical potential v — fi/ B. 
The integrals over are given by 

/C30 
dKy i/„(C - Ky) H^ie ~ Ky) e'^V i^^^ ) - ^ (i-yf - ^ ^ ^ (ij) 
-OO 

JtnAC. C, t) = / dK,e{v-\Kl-n- i) e^--iC-0-¥^lr (^8) 



The integration over the transverse wavevector can be carried out by using the identity 
/ dk Hnik + a) Hn{k + b) e-''^ =2''^/^n\ Ln{-2ab) (19) 

J OO 
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for arbitrary complex a and b. One finds: 



(20) 



with the two-dimensional vectors pj^ = {^,ri) and = (^',?y') The integral over the 
longitudinal wavevector can be rewritten in terms of Fresnel integrals [S] : 



1) Jl eiKC-Cr/r p 



2t 



•^||,«,i.(C,C',t) ^d{y-n 
with: 

F{a, b) ^ C{a + b)~ iS{a + b) + C{a - b) - iS{a - b) 
Combining lfTC|l with and (|^ . we get 



(21) 
(22) 



G^(r,r',i) = 



B3/2 



■KP-L-pl)' + 5i(C+«')('?-';')+ii(C-C')VT 



47^3/2^-1/2 

x^'e-'("+^)^i„[i(p^-pl)^] F 



y/T{v-n - i) 



\/2^ 



(23) 



The prime indicates that the summation is only over those values of n for which 
— (n + i) is non-negative, i.e. over the Landau levels that are at least partially 
filled. It should be noted that the right-hand side is continuous as v equals a half-odd 
integer, since F(a, b) vanishes for a = 0. 

The function G(r, r',i) can be evaluated along similar lines. Instead of 
>^||,n,j/(Ci C'l ''') 0116 encounters an integral like 118(1 . without the step function, which 
is easily evaluated. As a result we get 

G(v,v',t) = ^-'^ ~ i^'^T Q-j{p±-p'±?+¥(^+i')(n-v')+h'^(C-C'f/r 



4,^3/2^1/2 

oo 

X E e-'("+^)^i„[i(p^-pl)^] 



(24) 



n=0 



Here the summation extends over all non-negative integer n. We should add a small 
negative imaginary part to r to ensure convergence. A simpler form for G is obtained 
by using the identity [H| 



for |6| < 1. In this way (|24|) becomes 
G(r,r',i) 



(25) 



(1 + i)B3/2 l^)(p^_p,^)2_^li(^^^,)(^_^,)^l;(^_^,)2/^ 



87^3/2^1/2 sin(iT) 



(26) 

Clearly, G is singular for t equal to a multiple of 27r. 

Information on the equal-time correlations is obtained by putting t ~ Q. The 
function G is equal to 5{y — r') for < = 0, as is obvious from its definition. To 
determine the static limit of G^ one uses the asymptotic expansion ,8j for u — > oo: 



C(M)-iS'(u)«i(l-i) + 



'27r u 



(27) 
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In this way one finds from H23|) 

^3/2 



G^(r,r',i = 0) = 



-3(P_L-pl)' + 5i(?+C')(')-')') 



e 4 



x^'L„[i(p^-pl)2] sin[v^ 



^)IC-C'I 



(28) 



This result agrees with that found before |2j, if spin degeneracy is taken into account. 

The behavior of the time correlation functions depends on the orientation of the 
position difference vector with respect to the magnetic field. In particular, one may 
consider purely transverse and purely longitudinal directions separately. 

For a purely transverse position difference, e.g. for S, C with rj — rj' and C = d 
one finds from l|23(l : 

R3/2 „ , 

n 

X C (v/t(!/ - n - i)) - 15 (v/T(i^-n - i) 



(29) 

This expression gets particularly simple if only the lowest Landau level is filled. For 
I < < f one gets: 

^3/2 



G^(r,r',i) 



2^3/2 



C 



^r{.-\))-isUr{.-\) 



(30) 



The spatial decay is Gaussian, with a Laguerre-polynomial modulation in the general 
case H29|l . As a function of time, is algebraically damped, with a periodic phase 
factor and a Fresnel-integral modulation. For large t the latter drops out, as follows 
from (|27|l . Owing to the Laguerre polynomials in (|29|l . the asymptotic form of the 
coherent structure function S'co/i(r, r', i) for large t varies with the magnitude of the 
transverse position difference, even if the Gaussian prefactor is regarded as a trivial 
normalization factor. As we shall see, this feature disappears when the position 
difference becomes purely longitudinal. 

The behavior of G for transverse separations is obtained from H2()|l with rj = rj' 
and C = C inserted: 

i)B3/2 



G(r,r',t) = 



(1 

37^3/2^1/2 sin(i 



-e4 



,COt(iT)(J-^')^ 



(31) 



The function depends on the position difference through a phase factor only. As to 
the time dependence, the same singular behavior for r — 2'Km as in H26() is obtained. 

For a purely longitudinal position difference, the behavior of G^ is found from 

by taking ^ = ^' and rj — ij' . One gets: 



G^(r,r',i) = 



^3/2 



47^3/2^-1/2 



e2 



|i(C-C')Vr 



F 



c-c 

x/27 



(32) 



The spatial dependence is no longer Gaussian: it is contained in a phase factor and 
in the Fresnel integrals. As a function of time, a qualitatively similar behavior as in 
the transverse case is found, with an algebraic damping, a periodic phase factor and 
a Fresnel-integral modulation. For large t the latter disappears, as before, so that 
then Gf^ depends on the position difference through a phase factor only. Hence, for 
any value of the longitudinal position difference, Scoh{i^, r', t) has the same asymptotic 
form for large t. 
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Figure 1. Bulk coherent structure function for v = 2 (a) and f = 5 (b), as a 
function of time t. The curves, which are normalized to 1 for t = 0, represent 

the structure function for coinciding positions ( ), for a longitudinal position 

difference |C~C'I = 1 ( )i fo'' ^ transverse position difference IC — = 1 
( )• 



The behavior of G for longitudinal separations immediately follows from H26|l as 

As in the transverse case, G depends on the position difference through a phase factor 
only. The time dependence is dominated by the usual singularities. 

Finally, let us consider the time dependence of G^ and G for coinciding positions 
(r = r'). It is obtained by putting ^ = ^' in ((23) and ifH)) : 

G^(r, r, t) ^ ^'e-'("+i)- [g [^/riv-n-^)) iS (Vri^^n^^,) 

n 

(34) 

(1 + 1)5^/2 

G{r, r,t) — — . , . — - (35) 
STT-'/^r-^/^ sm(ir) 

For large t, the function reduces to 

n \ 2 y 

The sine function in the numerator, with rig the label of the highest occupied Landau 
level, ensures that is non-singular for r — 2Trm, in contrast to G. 

In figure n we show numerical results for the normalized bulk coherent structure 
function Scoh{^,^' ,t)/ Scoh{^,^' ,0), at two values of the reduced chemical potential i/. 
The anisotropy in the behavior of the structure function for non-coinciding positions is 
clearly visible. The curves that correspond to purely longitudinal position differences 
show characteristic revivals around r = 27rrn, with integer m. The revivals are less 
prominent for purely transverse position differences. These features can easily be 
understood from the interference effects in the sums in (|29|l and (|32|l . 
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4. Bulk intermediate scattering functions 

The intermediate scattering function in the bulk foUows by Fourier transform with 
respect to the spatial coordinates. The incoherent part is obtained by inserting (|^ 
and 1123 in ® with ®: 

5„,(k, t) = J d{p - p') c-'^^ iP-p') G;ir, r',t) G(r, r', t) 



327r''r sm( ir) ^-^ " 



(37) 



with K = Kz) = 'k/\fB and p = {pj_, () = ^/B r. The integrals are defined as 
^j_,n('«±, t) = J dp^ exp [ipi/(e"'^ - 1) - ik^ ■ p^] L„ {\pi) (38) 



J\lnA>^z.r) ^ JdC e-'^-'^F* 



0, 



c 



/2t 



(39) 



To evaluate J±^m we first perform the angular integral in terms of a Bcsscl function. 
Subsequently, we use the identity [7]: 



y 



2a(l - a) 



(40) 



with Rea > and y > 0. In this way we get 

Ji,„(/«i, r) = 47ri e-'("+3)" sin(ir) exp e"'") 4] L„ [2 sin2(iT)4] (41) 

For the integral over the parallel coordinate we employ the identity 



dz e 



-ikz 



F{a,zb) 



for real fc, a and 6, which follows by a partial integration. In this way we find 

2(1 + i) 



sm 



v/2(i/ - n - i)K^ 



(42) 



(43) 



Combining the above results, we have found the incoherent intermediate scattering 
function as 



27r-^KzT 

sin ^2(i/ - n - i)KzT L„ [2 sin^(iT)K5^] 



(44) 



For the special cases of purely transverse or purely parallel wavevectors the incoherent 
intermediate scattering function can be simplified. For = it becomes an 
undamped periodic function of r with period 27r, whereas for k_l = it is a decaying 
function as in the general case. 

For t = {] the scattering function reduces to 



^3/2 



S.„,(k, t = 0) = ^Y. V^i'^-n-k) 



(45) 



which is equal to the bulk particle density Hoo at the chosen chemical potential. The 
time derivative at i = satisfies the sum rule 



d 

-g^Sinci^, t) 



(46) 



t=o 
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We now turn to the coherent intermediate scattering function. By substituting 
in iP with we get 



ScohO^^t) = -i?-3/2 J d{p-p') e--(^-''') |G^(r,r',t)|' 



167r37 



The integrals are defined as 



^||,™„,.(«.,t)- J dC e--^?F* 



0, 



c 



/2r 



c 



/2t 



(47) 

(48) 
(49) 



The angular integral in J±^mn leads to a Bessel function as before. The radial integral 
can be determined by using the identity '7' for v > 0: 



„ '4k1 r-"i+"/l,^2\ 7--n+m/'l„2\ 



(50) 
(51) 



= (-ir 

which yields 

The two Laguerre polynomials are not independent, as they are related as |10| 

(-1)" m! 2-" L,7;"+"(z) = (-1)" n! z"" L;^"+™(z) (52) 

Hence, JL,mTi('^_L) is positive for all m, n > 0. For the integral over the parallel 
coordinate we use the integral identity (|A.1|I from the appendix. With its help one 
gets 

>^j|,mn,iv('*z: ''') ~ 

2i 



I 

2i _ii„2 r 



wl - + WmY 



Wn 



— C 



A) 

wl - (k^ - Wmf I (53) 



with the abbreviation Wj = ■\/2(j/ — j — i). 

Substituting and (|5^ in l|17|) and employing the symmetry with respect to 
an interchange of the summations over m and n, we obtain the coherent intermediate 
scattering function in the form 

ScoHiKt) = e-i«i (-1)"^" (I'^i) (^4) 

2 



27r2Krr 



|sin [{^kI - K^Wm + n - m) r] 6* wl - {k^ - w,„) 
- sin [(^K^ + K^Wnj + rt — m) r] 9 — (k^ + w^)^ | 



(54) 



For vanishing k± the product of the Laguerre polynomials equals Smn- As a 
consequence, the coherent intermediate scattering function gets the simpler form: 

^3/2 , 

S'coh(k, i) = — — j- sin [(i/i^ - t] 6'(2w„-|k^|) (55) 

ZiT \Kz\T 
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A second special case arises for = 0. Expanding the factor between curly brackets 
in (|54|l for small one finds: 

S^ouiKt) = e-^'^i Y^Y^ (-1)™+" L-r+" ihl) ihl) 

ni n 

X {cos[(m - n)T] [(1 - dmn) 0{m - n) + (1 - (5m„) 9{n - m) Wn] + S^n Wm] (56) 

In this case the coherent intermediate scattering function is an undamped periodic 
function of r. 

The static function follows by putting < = in (|54(l . Rearranging the terms one 
may rewrite it as 

R3/2 „ I I 

5,„,(k,t = 0) = e-i'^- E (-1)"^" (^4) i,T"+" (14) 

X [{Wm +Wn - \Kz\) (w,„ + W„ - |k^|) - (Wm - W„ - \Kz\) 9 {Wm - Wn - \Kz\) 

- {Wn - Wm - \Kz\) {w„ - Wm ~ \Kz\)] (57) 

In contrast to Sine, the coherent intermediate scattering function is an even function 
of t, as is clear from (I54II . Hence, the time derivative of Scoh at t = vanishes, whereas 
the time derivative of Sine at t = is different from 0, as we have seen in H46|l . 

On comparing H44() and H54() one notes that the two parts of the intermediate 
scattering function are given by rather different expressions. However, they can be 
made more analogous by employing the identity 

C30 

(-1) e ^ L,„ ^ (2K_Lj i« (2^-Lj 

m=0 

= (-1)" e-'("+5)^ L„ [2sin2 (It) k^] exp (ie"'" kI) (58) 

which may be proved by using H5U|) in reverse order, and then H25|) and H4U|) . In this 
way we may write (|44ll as 



oo 



m— n 

xir„™+" (i^i) sin(«;„«:,T) (59) 

This result could have been obtained in an alternative way by substituting H24|l 
(instead of H2t)|) ') in ^ with and taking similar steps as in (|47|I - H54|I . 

5. Bulk dynamical structure factor 

The dynamical structure factor follows from the intermediate scattering function by 
taking a Fourier transform with respect t. For the incoherent part of the structure 
factor we find from l|59|l : 



/oo 
dT e'^- S,nc{k, 
-oo 



t) 



dI/2 ^ 

e"^''^ E E'(-i)'"+" (^4) L--^"" ih^l) W - ^muH] (60) 



' ' m— n 



where we introduced the rescaled frequency zu ^ u/ B , and the abbreviation 

Vmnit^) = n + i + ^K~^ {^kI - w + m - (61) 
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The dependence of Sine on lo is determined by the step function. For negative lo it 
imphes that the summation variable m is confined to values below v — as is n. For 
positive LO this need not be the case. 

The coherent part of the dynamic structure factor follows from H54|l as 

5..(k,.) = -|;^e-^''iE'E'(-ir 



xL-'"+" (14) L,;"+™ (14) e {v z.„,„(-|c.|)] (62) 

Again the Ci;-dependence shows up through the step function only. 

From (|52ll it follows that the incoherent part of the structure function is positive 
for all k and cj, whereas the coherent part is negative. Furthermore, on comparing (|6U|I 
and H62() one finds that the two parts of the dynamic structure factor are connected as 
S'co?i(k, w) = — S'mc(k, — |a;|). This relation is in agreement with the general identity 
H12|l . since the structure factor is symmetric under spatial inversion. It implies that 
the total dynamic structure factor at T = vanishes for negative io. 

If the wavevector is purely transverse, the expressions H6U|I and H62|l cannot be 
used as such. It is more convenient in this case to take the limit = in ()59|l and 
to use (|56|l before taking the Fourier transform. In this way one gets 

Rl/2 °° I 

7n—0 11 

xL-™+" (14) (14) u)„ 5{vj ~m + n) (63) 

and S'co/i(k, w) = — 5mc(k, — |w|), as before. Both parts of the dynamical structure 
factor are sums over delta functions, with peaks at the integer resonance frequencies 
m — n. 

For wavevectors parallel to the field one finds from H60() : 

^1/2 , 

5„c(k,C^) = V e[v-Vnn{u:)] (64) 

' ' n 

and an analogous expression for the coherent part. The sum of step functions reaches 
its maximum value for w = ^k^.. 

The numerical results in figure[21show how the coherent dynamic structure factor 
changes from a function with steps for a purely longitudinal wavevector to a set of 
discrete spikes for a (nearly) transverse wavevector. A similar behavior is found for 
the incoherent dynamic structure factor, as shown in figure |3| For this case the curves 
are asymmetric with respect to a change of sign of ui. 

The sum rules (145(1 and ((46|l for the incoherent intermediate scattering function 
yield 

/>oo />oo 

du) uiSincO^,!^) = Trfc^rioo (65) 

The above results for the bulk dynamical structure factor at T = can easily be 
generalized to finite temperatures. In fact, for the incoherent part (|ll|l . one may use 
an identity analogous to (jSJ to obtain the expression for finite T from that for T = 0. 
Since the chemical potential in ((60|l only occurs in the step function, the integral over 
jj,' is trivial. In this way one easily arrives at the result for finite temperature: 

I I m=0 n=0 




Figure 2. Bulk coherent structure factor Scoh(^y'^) / B^^'^ for u = 5 and |fi| = 1, 
as a function of ro. Tlie curves represent tlie structure factor for 5 = (a), 6 = 7r/4 
(b), and 9 = 1.55 (c), with 9 = arctan(|K.^ 




-4 -2 2 4 -4 -2 □ 2 4 -4 -2 2 4 



Figure 3. Bulk incoherent structure factor SincO^, lo)/B^/'^ for v = b and |k.| = 1, 
as a function of ro. The curves represent the structure factor for 5 = (a), 9 = tt/A 
(b), and 9 = 1.55 (c), with 9 = arctan(|K^|/K2). 



xL 



— m+n / 1 _2 



(i-i) L--^- ih-l) ^ 1 (66) 

where the temperature is given in units of B. The coherent part of the dynamic 
structure factor for finite temperatures foUows from H12|l : 



oo oo 

1 .-2 



R1/2 

C n, , .\ \^ T\m+n T -m+n r 1 2 \ t -n+m f 1 2 \ 



' ^' m=On=0 



11 11 



(67) 



6. Vanishing magnetic field 

In the Umit of vanishing magnetic field the resuhs of the previous sections are expected 
to reduce to the well-known expressions for a non-magnetized degenerate free-electron 
gas. We may check this by considering the limiting behavior of the functions en G 
in position space. 

The expression (123(1 for contains a sum that extends to infinity for _B — s- 0. 
Hence, the order n of the Laguerre polynomial may become large. On the other hand, 
its argument gets small, since the position variables scale with ^/~B. In this regime the 
Laguerre polynomials may be approximated as 

L„(u) « e^" Jo (V2(2n+ 1)m) (68) 
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Putting nB — v, we may write the sum in (|23|l as an integral for small B. Replacing 
the scaled variables by the original ones, we get for B = 0: 



G,{r, r', t) = ^^5^ e^'(^-^')V* J^' dv e"-* Jo (%/2^ |rj_ - 
With the use of the relation (jA.7|l of the appendix we obtain 



(69) 



for B = 0. For large t the second term dominates. It factorizes in a space-dependent 
and a time-dependent factor. The latter is proportional to t~^ (apart from a t- 
dependent phase factor). This is in contrast to the behavior in the case with field. 
For r = r' we have seen in H36II that Gfj, has a decay proportional to t^^^^ times the 
ratio of two time-dependent sine functions. 

The limiting form of G for vanishing B follows immediately from l|26|l. One finds 

G(r,r',t)^-j-l±Lje^'l-'lV* (n) 

It decays algebraically in time, and is free from the singularities that characterize the 
expression (I26|l for the case with field. For small fields the poles in H2()l) at r = 27rm, 
or < = 27rm/_B, with integer m ^ 0, shift towards oo, and disappear for vanishing field. 
Only the singularity ai t = remains. It is connected to the delta function singularity 
of the static function G(r, r', t = 0). 

The incoherent and coherent parts of the dynamic structure factor follow by 
taking Fourier transforms, as we have seen in the previous sections. As a useful check 
of these results we may verify that the limits of H60|) and H62|) for vanishing field agree 
with the well-known results for a field-free degenerate fermion gas. Starting with the 
incoherent part we first use H5()(l in reverse order. The resulting integral over u, 
with V — contains a Bessel function Jq(uk±). For small field, k± becomes large so 
that effectively only small values of u contribute in the integral. Hence, one may use 
H68|l for the Laguerre polynomials in the integral over u, so that the latter becomes an 
integral over the product of three Bessel functions. It can be evaluated as a product 
of two inverse square roots ,7^- In this way we have found for large v, and large m, n: 

« ^ e [Wmn - {v^ - 2m - 2n f] (72) 

TT^yiGmn — (ti2 — 2m — 2n)^ 

Employing this asymptotic relation in H60(l and replacing the sums over m, n with 

integrals, by writing x — mB and y — nB, we get for B = 

I r-=° rt^ 1 
S'j„c(k, w) = — r / dx dy 



x( 



2^2|fc,| Jo Jo ^-x^ + 2xy -y^ + kl{x + y)~ ifc^ 
[-x^ + 2xy -y^ + kl{x + y)~ ^ki] 

„2 I o^„, „,2 I „,\ I n, „,\ ('11,2 , ,\2 



xe* -2;^ + 2xy - y^ - ki{x + y) + 2uj{x ~ y) - {\ki - loY + 2^ikl\ (73) 
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The two step functions determine a region in the xy-plane that is enclosed between 
two parabolas. Integrating first over x -\- y, and subsequently over x — y, we arrive at 
the expression 

2" 

, / Ll) An \ 

1 



27rfc 



(74) 



valid for —\f2\l k + ^k^ < uj < y/2jl k + ifc^. The coherent part of the dynamic 
structure factor may be discussed along similar lines. One finds in the field-free limit 



S„»(k,..) = -^ 



(75) 



for k < 2-^2/7 and \uj\ < y/2jl k — ifc^. The expressions (|74|l and (|75|l agree with the 
well-known results for a degenerate fermion gas without a magnetic field [3], [Hj, 

7. Edge effects: hard wall perpendicular to the field 

In the presence of a hard wall at z = the energy eigenfunctions ipn ,nj,n(r) in the 
domain z > are related to (|14|l by a refiection principle: 

(fny^n^A'^) = 4s ['0n«,n.,n(r) " "^n^ ,n (? )] (76) 



V2 



with f = (x,y,—z). Hence, both and G can be obtained immediately from the 
results of section 3. We will only consider the case z = z' . From H23|l we get: 

G,(r,r',i) = e-i(''---^)^ + ^HC4-«')(.-.')^VH«+^)^ L„ [i(p^ -pl)^] 



X <^ F 



VT(i^-n - i),0 



(77) 



G(r,r',i) 



^^/2 gji^l' 1^) V / 



Likewise, we find from H26() : 

'87r3/2Ti/2 sin(iT) 

For equal position vectors (r = r') the Laguerre polynomial and the space-dependent 
exponential drop out from (|77|l . so that one gets: 

^3/2 



G^(r,r,i) 



47^3/2^-1/ 



X <^ F 



1/2 



v/T(zy-n - i),0 



e2iCVr^ 



v/T(i/-n - i), ^/-C 



The expression (|78|l for G reduces in this case to 

(l-Ki)B3/2 



G(r,r,t) 



87,3/2^1/2 sin(iT) 



(79) 



(80) 
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G^(r,r',i = 0) 



^3/2 



27r2 



e 4 



sm 



2y2( 



ly — n ■ 



(81) 



This expression may also be obtained directly from H28|) . For coinciding positions the 
Laguerre polynomial and the exponential prefactor again drop out. The function G 
reduces to d{r — r') in the static limit, as in the bulk case. 

The expressions (I77|l and H78|) (or H79|l and H8UI) for coinciding positions) for 
and G may be compared to their counterparts and (or and (|^ ) 
for the bulk. The contributions from the reflection in the wall lead to a change in 
the dependence of the correlation function on time. This becomes manifest upon 
choosing positions near the wall. A Taylor expansion for small C leads to the following 
expressions for the correlation functions near the wall: 



:/-2 r3/2 

G,ir,r',t) « - i,^,. e-^^"--''-) +^'«+«')(''-''') 



7^3/2^3/2 



and 

G(r,r',i) 



-iS i^\/T{iy - n - i) 

(1 -i)C2B3/2 



- \l - c 

TT 



Vr(z^-n - i) 



(82) 



47r3/2^3/2 sin(ir) 



iicot(ir)(p_L-pi)' + 5i(«+«')(f)-V) 



(83) 



For equal positions these become: 

1(2^3/2 



G^(r,r,i) 



and 



^3/2^3/2 

n 

-i^(v/^ 
G(r,r,t)« 



^' {e-("+^)- [g (Vr(.-r.-i)) 



(1_ 1)^2^3/2 



(84) 



(85) 



47r3/2T-3/2 sin(ir) 

For large r the Fresnel integrals between the curly brackets in H82|l and H84|l can be 
neglected, so that both expressions show an algebraic decay proportional to r"^. From 
H82|l we get in this way 



G^(r,r',i) w ^^^^^^^ Q-k(.P±-p'±f+Hi+i')iv-v')-i'^T 



c^'L„[i(p^-pl)2] V2(i^-n-i) 



(86) 
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Figure 4. Coherent structure function as a function of time t, at !^ = 5, for 
position difference |^ — = 1 (a), and for coinciding positions (b). The curves, 
which start at 1 for r = 0, give the normalized structure function for distances ^ 
from the wall equal to 0.1 ( ), 1.5 ( ) and oo (bulk, — • — ). 



which for coinciding positions reduces to: 

(r, r, t) « e-"-^ V ^2{u-n--^) = ^ e""^^ (87) 

n 

with Uao the bulk particle density. 

Comparing H87|) with (|36|l . one concludes that the behavior of the correlation 
function G^, and hence of the coherent structure function Scoh, has changed owing to 
the presence of the wall. In fact, whereas in the bulk the coherent structure function 
for coinciding positions decays proportional to l/r for large time (with a periodic 
modulation), the decay near the wall is faster, namely proportional to r~^. The 
function G has changed as well, with a decay proportional to t^"^/^ near the wall, 
which leads to a corresponding change in the behavior of Sine- However, the recurring 
singularities at r = 27rm dominate the picture completely in this case, and these are 
not altered by the wall. 

In figure0]the normalized coherent structure function S'coh(r, r', t)/S'co/i(i", r', 0) 
is shown for several values of the distance from the wall C, and the position difference 
1^ — ^'|. The faster time decay of the structure function for positions near the wall is 
manifest. Revivals are (almost) completely suppressed in this case. 



8. Edge effects: hard wall parallel to the field 

Let us now consider wall effects in a different geometric arrangement, namely for a wall 
parallel to the field. For a hard wall at a; = the energy eigenf unctions '4}ny,n^,n{Y) 
for X > can be expressed in parabolic cylinder functions [H] as 



{r^c ^?„K)-i/2 [V2(^'-^^.)]} 



As for the bulk case, the wavevectors are given as = 27rni/(Lvi3), with integer ni 
for i — y,z. Furthermore, the function £„(/«) occurring in the index of the cylinder 
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functions is determined by [2] 

^e„W-l/2(-V2K) =0 (89) 

The same function shows up in the eigenvalue associated to H88|l 

En^,n^,n = Ben{Ky) + ^kl (90) 

The time correlation function for the degenerate case follows by inserting (|88|l 
in 0. We will only consider the case x = x' . Using (|18(l and (|21(l . with n + ^ replaced 
by £„(Ky), we get 

o3/2 ^ , pea 

(1 (r r' t\ — p^i(C-C') /r \^ / iKj,(r;-77')-ie„(Ka)r 

P { ~ £n{Ky)], ] (91) 



where Kn{v) is defined by the condition e„[K„(i^)] = with n < i^ — i. In an analogous 
way, the function G for a; = x' is found as 

(l-i)S3/2 



iKy(r;— r;') — ie„(/ty)r 



^el(.„)-l/2[^/2(C-.,)] 

r^e"^LK)-i/2[V2(e'-«,)] 

Here, the sum extends over all non- negative integer n and the integral over all Kj,, in 
contrast to the sum and integral in H91|l . 

For large values of x the expressions (|91|l and (|92l) reduce to their bulk forms (|23|l 
and H24I) . In fact, for ^—+00 only large values of Ky contribute to the integrals. For 
these large Hy the function en{ny) reduces to n+^. Furthermore, the parabolic cylinder 
function £'„(\/2it) for integer n can be rewritten in terms of a Hermite polynomial 
as 2~"/^e~" /^Hn{u). Finally, the lower bounds of both the integrals over ^" in the 
denominators in (|91(l . (|92|l and of the integral over Ky in (|91|l can be extended to —00. 
Evaluating the integrals with the use of (|19|l . one recovers H23II and (|24|l . 

As a special case we consider the time correlation functions and G for equal 
position vectors, as in section 3 and 7. One finds from 19111 and (|92|l : 



, poo 

G.ir,r,t)^^^^^jj^^ J d.ye 



|g ^y'r[i^-£„(Ky)]^ - 15 ^y'T[j/-£„(Ky)] 



(93) 



and 



^el(.„)-l/2[^^(g-^.)] 

r^e'i?L(«j-i/2[V2(r--.)] ^ ' 

These expressions are the counterparts of (|34f) . H35|l for the bulk, and of (|79|l . H8U|) for 
the other geometry. For large t the Fresnel integrals in (|93|l can be replaced by their 
asymptotic value i. 
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Figure 5. Coherent structure function as a function of time t, at v = 5, 
for position difference \p — p'\ = 1 (a-c), and for coinciding positions (d). 
The first three correspond to 9 equal to (a), 7r/4 (b), 7r/2 (c), with 9 = 
arctan(|?7— ?7'|/|f — f'l). The curves, which start at 1 for r = 0, give the normalized 

structure function for distances ^ from the wall equal to 0.1 ( ), 1.5 { ) 

and oo (bulk, — ■ — ). 



In figure[3nunierical results for the normalized coherent structure function, which 
follows by inserting (|91|1 or (|93|l in are presented. Curves are drawn for several 
values of the distance ^ from the wall, and for various position differences, with a few 
different orientations of the latter. The general picture of these figures is the same as 
that obtained for the other geometry: near the wall the time decay of the structure 
function is considerably faster than in the bulk. More details on these figures will be 
given below, when discussing the time correlation functions near the wall. 

For t — Q the correlation function can be simplified by using the asymptotic 
relation (|77|l . One gets from 

R3/2 , poo 

G,(r, r', t = Q)^ I 



^ ' ^ (95) 



which agrees with a result from j2j, if spin degeneracy is taken into account. As before, 
the second correlation function G(r, r', t) becomes a spatial delta function in the static 
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limit. 

As in the previous section, we want to study the behavior of the time correlation 
functions in the vicinity of the wall. To that end we start from ()91f) and (|92|) . and use 
the identity 



dn 



r<'^??„(.)-i/2[V2(e"-«)] 2vr 

(96) 

Furthermore, we make a Taylor expansion of the parabolic cylinder function in the 
numerators of 1911) and H92|l . Subsequently, we employ the Wronskian identity [S] 



r[-£„(«;) + i] 



(97) 



which holds as a consequence of H89|l . Upon changing the integration variable from 
Ky to e in H91II . we finally obtain the time correlation function for positions near the 
wall as 



G^(r,r',i) 



f2 r3/2 



27r3/2^1/2 



c-c 

V2^ 



Likewise, (I92II gives 

(1 - ^)eB'^' ii(c-c')v. 



G{r,r',t) 



27r3/2-ri/2 



e2 



— ri — ler 



(98) 



(99) 



,1=0 "'"+1 

These expressions may be compared to (|82|l and (|83|) for the other geometry. We 
will now consider the time correlation functions in the vicinity of the wall for some 
particular orientations of the position difference, and for the case of coinciding 
positions. 



8.1. Longitudinal position difference (y = y', z ^ z' ) 

Putting r] — rj' in the expression l|98|l for , we can evaluate the integral by a partial 
integration. The result is 



c-c 



sm 



For large separations |C — C'l we obtain with the use of (|27|> 

^2 ^3/2 



(100) 



G4r,r',t) 



)cos ^2{v-n-^) IC-C'I (101) 



TT^IC-C'P 

For i = the expressions (I100|l and p01|) reduce to static results found before [2- For 
large r the second term between the curly brackets in (|100|) dominates, so that one 
gets: 



;f2 o3/2 , 

G,(r,r',t)«^^^e-'-5] 



^2(1^-71-1) IC-C'I 



IC-C'I 



(102) 
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Upon comparing this result with (|86|l , which holds for the other geometry, one finds a 
similar behavior as a function of time: both have a decay proportional to t~^, which 
differs from the t~^/^ decay of for the bulk. These changes in the behavior of 
entail modifications in the decay properties of the coherent structure function that 
have been discussed in the previous section. The differences in the spatial dependence 
between (|86|l and (|102|) are analogous to those found in the bulk static correlation 
function H28II for transverse and longitudinal directions of the position difference, 
respectively. 

Numerical results for the coherent structure function with a purely longitudinal 
position difference are given in part (a) of figure |5l For a small distance from the wall 
= 0.1) the structure function decays rapidly. It does not show the revivals around 
T = 27rm, which are a prominent feature of the bulk structure function, and which we 
have seen already in figure ^ For intermediate distances = 1.5 in the figure) these 
revivals have lost part of their strength. 

As to G, the integral in (|99|l becomes trivial for rj = rj' , so that one gets 

^ ' ' ^ 4^3/2^3/2 sin(iT) ^ ^ 

which may be compared to (|83|l for the other geometry, and to H33|l for the bulk. As 
before, the time dependence shows the usual characteristic recurring singularity and 
a faster decay than in the bulk. 

8.2. Transverse position difference (z = z', y ^ y' ) 
In this case, the expression (|98|) for G^, becomes 

c2 o3/2 , i-y 

G.ir, r', t) « de e'-(^)(''- 

X [C yT{v-e)\ - xS yr{v^e)\ } (104) 

In general, the integral cannot be simplified further. However, for large |?7 — r/| one 
may find its asymptotic behavior, as it is determined by the upper boundary. One 
gets 



e2 



^3/2 gfirisgn (»)-r/) , 



V27r3/2 |?7-77'|3/2 ^ 



dv 



-3/2 



(105) 



which for t = reduces to an earlier result 2 . As in p01|l the time dependence is 
given by a trivial phase factor. The space dependence is quite different from that of 
itniTll . 

For G as given by H99|l . the choice C = C' produces only a slight simplification. 
The integral in H99|l remains as it stands; it is difficult to evaluate, even for large 

In part (c) of figure [S] some numerical results for the coherent structure function 
with a purely transverse position difference are presented. As before, the structure 
function decays rapidly, if the distance from the wall is small (^ = 0.1). The behavior is 
rather similar to that for the other geometry (see part (a) of figureEj). At intermediate 
distances (^ = 1.5) some revivals are present, as in the longitudinal case (see part (a) 
of figure El, albeit with a smaller amplitude (note the difference in scale). At this 
distance from the wall the structure function is quite different for the two geometries, 
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as can be seen by comparing with the corresponding curve in part (a) of figure^ In 
the bulk the revivals have a more complicated structure, as we have seen already in 
figure ^ and in figure^ The intermediate case of a position difference that is oriented 
at an angle 7r/4 with respect to the field leads to a behavior of the coherent structure 
function that interpolates between the purely longitudinal and the purely transverse 
cases, as can be seen in part (b) of figure El 



8.3. Coinciding positions (y — y' , z = z' ) 
For coinciding positions we get from l|l(K)|l 

1^2^3/2 



G,(r,r,t)«-i|^5:'{e-("+^)^ [c {^r[. - n - ^)) 

n 

-iS (Vr(^^-n-i))] - e-'"^^ - - |) | (106) 



Remarkably enough, this expression is identical to H84I) . apart from the interchange 
of ^ and C- Hence, near a wall the time correlation function for equal positions has a 
universal form, independent of the orientation of the wall with respect to the field. 
The same conclusion is reached for G. Indeed, from l| 103(1 we find for ( = d^': 

^^-'-'^)-- j3/V3/2sin(^.) (10^) 

which has the same form as ()85|l. 

For large t the second term between the curly brackets in H106|) dominates. 
Introducing the bulk particle density, we may write for coinciding positions near 
the wall in this asymptotic time regime in a form analogous to (|87f) : 

G,,(r,r,i)«^e-''^^noo (108) 

T 

In accordance with the more general result discussed below (|106f) . the long-time 
asymptotic behavior of the equal-position time correlation function near a wall is 
independent of the orientation of the wall with respect to the field. The asymptotic 
form has an even stronger universality property: it is independent of the strength of 
the field as well. This follows by reintroducing the original (non-scaled) space and 
time variables, and the chemical potential /i. In terms of these variables (|108|l reads: 

2iT2 

G^(r,r,i)« — e-'^*no, (109) 

The right-hand side is indeed independent of the magnetic field (apart from the implicit 
dependence through the bulk particle density). It should be noted that in the bulk 
the long-time tail in G^ for equal positions does depend on the magnetic field, as we 
have seen in 

The universality property found for G^ does not hold for G. In fact, reintroducing 
the original variables in H107|l we get 

Clearly, this result depends on the magnetic field strength. It shows the same 
characteristic recurrent pole structure as we have found in the bulk (see H35|) '). As 
in the bulk, the poles at < = 2'iTn/B move when the field strength changes. 
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Numerical results for the coherent structure function with coinciding positions 
are given in part (d) of figure [3 In the vicinity of the wall = 0.1) the structure 
function is (almost) identical to that near a wall perpendicular to the field (see the 
curve for ^ = 0.1 in part (b) of figure^l), in accordance with the analytical results (|84|l 
and (|106|l . For intermediate distances from the wall the curves for the two geometries 
are rather different, while they agree again in the bulk, of course. 

9. Concluding remarks 

In this paper we have studied time-dependent pair correlations in a quantum gas of 
charged free fermions in a uniform external field. These correlations can be described 
in terms of two correlation functions G'p_T(r, r', f) and G'(r, r',i), which together 
determine both the coherent and the incoherent part of the time-dependent structure 
function S'(r,r',t), according to (PJ-lEIl- Our main interest has been focused on the 
behavior of these correlation functions for the completely degenerate quantum gas at 
temperature T = 0. 

Owing to the presence of the magnetic field, the particle motion is quantized in 
terms of Landau levels with an associated cyclotron frequency. One of the aims of the 
present paper has been to show how the correlation functions are influenced by this 
cyclotron movement. The results and express the functions G^,t=o and G as 
sums over Landau levels, involving Laguerre polynomials for the spatial dependence 
and Fresnel integrals for the time dependence. Their time behavior is sensitive to 
the orientation of the position difference with respect to the magnetic field, as is 
demonstrated in (|29|I - H33(I . The time behavior of the ensuing coherent part of the 
structure function, as shown in figure^ is characterized by recurrent 'revivals', with a 
markedly increased strength of the correlation effects at values of the time t equal to 
a multiple of the inverse cyclotron frequency. These revivals are particularly clear for 
position differences in the longitudinal direction. The incoherent part of the structure 
factor is influenced even more strongly: it is singular for the values of t mentioned 
above, as can be seen from H26II . 

The time dependence of the correlation functions and the structure function can 
be translated to a frequency dependence. Passing via the intermediate scattering 
function S'(k, t), for which we have derived the expressions (|44|l and H54|) . we arrive at 
the dynamical structure factor S'(k, w), as given by (|60|) and Ht)2|l . The curves in flgure 
[3 and 121 show how both the coherent and the incoherent part of S'(k, w) get a spike 
structure for purely transverse orientations of the wavevector. The gradual change of 
these curves as the wavevector gets a longitudinal component can be inferred from the 
flgures. For a purely longitudinal wavevector a step structure is found, in accordance 
with 

The above results pertain to the bulk of magnetized quantum gases. As is well 
known, these systems show additional interesting physical phenomena near the walls 
of the containers to which they are conflned. A second goal of the present paper has 
been to investigate how the time-dependent correlations change in the vicinity of a 
conflning wall. In general, these edge effects depend on the orientation of the wall 
with respect to the fleld. They have been determined for two speciflc orientations, 
namely for a wall perpendicular to the fleld and for a wall parallel to the fleld. In 
the former case the cylinder symmetry of the fleld is not changed by the wall. As a 
consequence, the time correlation functions are easily obtained from those of the bulk, 
with results given in l(77jl - ((7H| . The curves of figure^ show that the main influence 



Time correlations in a confined magnetized free-electron gas 



23 



of the wall is a faster time decay of the correlations. For a wall parallel to the field it 
is a lot more cumbersome to determine the correlation functions, as the symmetry of 
the problem is now lost completely. The general results for the correlation functions 
are given in H91|) - (|92|l . The curves for the coherent structure function, which are 
presented in figure [SJ show that once again the correlations die out faster than in the 
bulk. In particular, the characteristic revivals at multiples of the inverse cyclotron 
frequency, which are a conspicuous feature of the bulk coherent structure function 
for longitudinal position differences, are almost absent for positions near the wall. It 
appears that the wall perturbs the regular patterns of the bulk motion so strongly 
that the recurrences associated to the cyclotron frequency do no longer influence the 
correlation functions. 



Appendix 



In this appendix we consider several integral relations that are employed in the main 
text. The first identity, which has been used in section 4, is 



dz e 



2i 



-\kz 



F*{a,zc) F{b, zc) 



" k 
2i 



{e''^"/^ e [^a'c^ -{k + 2hcf 



-ikb/c 



+ 1 g-iife^/c^ |e-ifca/c g ^4^2^2 _ ^ 2ac)2] - (}^-/- 9 [Ab^c^ -{k- 2ac)2] } (A.l) 

with the function F as defined in H22|) and with positive a, b, c and real fc 7^ 0. To 
prove it, one starts by using a partial integration to rewrite the left-hand side as 

i%/2c , |Jgi(a+.c)^ _ gi(a-.c) = j 



dz e 



-i{b-zcy 



'(a,zc)| 



(A.2) 



To determine the integral at the right-hand side, we have to evaluate integrals of the 
form: 

/•oo 



G{k,a,b) 



dz e' 



F{a, zb) 



(A.3) 



for positive a, b and real k. We differentiate this function with respect to k. By adding 
a suitable multiple of G itself we arrive at: 



-2b 



2 dG(k, a, b) 
~dk 



— ik G{k, a, b) 



^^-ikz+ib^z^^^ F{a,zb){kA) 



Integrating by parts and evaluating the ensuing integral the right-hand side becomes 
a combination of delta functions: 

- 2^/^ b e"'"' [5{k + 2ab) - S{k - 2ab)] (A.5) 

By combining (jA.4p and (|A.5p we arrive at a differential equation, which can easily 
be solved as: 



G{k,a,b) 



2tt 



c 4 



{Aa%^~e)+G{a,b) e-3 



(A.6) 



with an as yet unknown function G{a,b). To determine it we put fc = 0, and 
differentiate both sides with respect to a. Using the definition HA.3|I we find, 
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by differentiating F{a, zb) and performing the resulting integral, that G'(0, a, h) is 
independent of a for all positive a and h. As a consequence G(a, b) is independent of 
a as well. Finally, by considering the limit a — > oo (for /c = as before) in (jA.Sp and 
employing (|27|) . we prove limQ_+oo G(0, a, 6) = \/27r/&, and hence lima^oo (5(0, 5) = 0. 
From these results we conclude that the term with G in (jA.6p can be omitted, so that 
G{a,b) is determined. After these preparations it is an easy matter to demonstrate 
the vahdity of l|A.l|l . 

A second identity (used in section 6) is 



Jo 



dv e 



Jo iV^b) F (,/{^t,^] = -i e^'^^/* ^ ( V^, ^ 



93/2: 



'It 



\/a(62 + c2) (A.7) 



for non-negative 6, c (with VP + c? > 0) and positive a, t. Before starting with the 
proof of this identity, we remark that its dependence on a is trivial, as follows from 
a rescaling of the variables. Hence, we will put a = 1 in the following. The function 
F in the integrand at the left-hand side of \K.1\ may be replaced by the integral 
representation 

F ( V(l-«)t, - a/^c-^'^'/* dk 6--'=-"'=' (A.8) 

V 2VV V TT J-^T^ 

as may be established from the definition of the Fresnel integrals (or by comparing (|18|l 
and H21|l ). Inserting the integral representation in ljA.71) . we arrive at a double integral 
over V and k. Introducing new integration variables p and 9 by writing v — siv? 9 
and k — pcos9, we rewrite this double integral as 

2/" dp p^ e~'*P'' [ d9 sm9 Ja{bpsm9) e-'^P""^" (A.9) 
Jo Jo 

The integral over 9 turns out to depend on b and c only through the combination 
Vp + (?. This can be seen by using the standard representation for the Bessel function 
and writing the (^-integral as 

1 r'^ fiiT 

— d9 sin6»e-'^''™''™^'^-''='''=°''^ (A.IO) 

27r Jo Jo 

The argument of the exponential function is proportional to the scalar product of a 
unit vector in the direction of (0, tp) and a constant vector with cartesian components 
(6, 0, c). Since the integration extends over all directions of the unit vector, the integral 
l)A.10|l , and hence l|A.9|l as well, indeed depends on b and c through the length VP + c? 
only. Once this is known, we may evaluate (|A.9I) by choosing 6 = and replacing c 



by ^/W+<? in the answer. Since the Bessel function equals 1 for 6 = 0, the integral 
l|A.9|l can easily be evaluated in terms of Fresnel integrals, with the result: 



2i 

— e sm c - 



(A.ll) 



ct V2i3/2 ^2VtJ 

Generalizing this answer to the case 6 7^ in the way indicated, and substituting it in 
the left-hand side of (|A.7p . we arrive at the desired result. 
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